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We investigate the Jeans instability of a galactic disk embedded in a dynamically responsive dark halo. It is shown that the disk-halo system 
becomes nominally Jeans unstable. On small scales the instability is suppressed, if the Toomre stability index <2t is higher than a certain 
threshold, but on large scales the Jeans instability sets invariably in. However, using a simple self-consistent disk-halo model it is demonstrated 
that this occurs on scales which are much larger than the system so that this is indeed only a nominal effect. From a practical point of view the 
Jeans instability of galactic disks is not affected by a live dark halo. 
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1. Introduction 

The classical paradigm is that at given total mass of the system 
a galactic disk is stabilized against local dynamical Jeans insta- 
bilities, if it is embedded in a dark halo. This can be seen for 
instance from Toomre's Q T stability index for a stellar disk, 



Qi 



KCT d 



3.36 GS d 



(1) 



where k denotes the epicyclic frequency of the orbits of the 
stars, <r d is the radial velocity dispersion of the stars, and 
the surface density of the disk. G denotes the constant of grav- 
itation. If all other parameters are kept constant, but is low- 
ered, Qi rises and the disk becomes more stable against Jeans 
instabilities. The physical reasoning is that the self-gravity of 
the disk, which has a destabilizing effect, is reduced by the sur- 
rounding halo. Similarly the onset of non-axisymmetric coher- 
ent large-scale instabilities of the entire disk such as the bar 
instability was thought to be damped by a surrounding halo. 
Ostriker & Peebles (119731) showed in their classical numeri- 
cal simulations of the dynamical evolution of a self-gravitating 
disk that the bar instability could be suppressed, if the disk 
was embedded in a halo potential. However, modern high- 
resolution simulations in which the surrounding halo is treated 
as a dynamically responsive system, have shown that actually 
the opposite is true. Athanassoula (2002) showed that in her 
simulations of the bar instability the bar grows stronger, if the 
disk is embedded in a live dark halo rather than in a static 
halo potential. This was explained there and, particularly, in 
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Athanassoula (2003) as due to the effect of the live halo on the 
angular momentum exchange within the galaxy. First doubts 
about an entirely passive role of the halo were already raised 
by Toomre ( 119771 ). These findings were supported by theoreti- 
cal studies of the swing amplification of shearing spiral density 
waves which is also enhanced, if the disk is embedded in a live 
dark halo instead of a static potential (Fuchs 2004, Fuchs & 
Athanassoula 2005 ). 

In this note we return to the severe local Jeans instability of 
a self-gravitating disk and investigate the effect of the presence 
of a live dark halo. In the next section we demonstrate how 
Toomre's concept of the Qi parameter can be generalized in 
order to take into account the effect of such a halo. In the final 
section we discuss implications for realistic disk-halo systems. 



2. Modification of the Q T stability index 

We study the Jeans instability of an infinitesimally thin galactic 
disk using the model of a patch of the galactic disk developed 
by Toomre (fT964l , Goldreich & Lynden-Bell ( [1965] ) and Julian 
& Toomre (119661 ) (cf. also Fuchs 2001). The patch is assumed 
to rotate around the galactic center and the differential rotation 
of the stars is approximated as a linear shear flow. The surface 
density is assumed to be constant over the patch. Polar coordi- 
nates are approximated by pseudo Cartesian coordinates (x, y) 
with x pointing in the radial direction and y in the direction 
of rotation, respectively. Toomre (1964 ) has calculated the dy- 
namical response of the disk to a small 'ring-like' perturbation 
of the gravitational disk potential of the form 



<t>kexp/(«f + kx) 



(2) 
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by solving the linearized Boltzmann equation. The induced 
density perturbation can be written in the limit oj — > as 
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where Io denotes the modified Bessel function and the back- 
ground surface density of the disk, respectively. In deriving 
eq. (O a Gaussian velocity distribution of the stars with a ve- 
locity dispersion cr d has been adopted. The disk is assumed to 
be self-gravitating, so that the density - potential pair has to 
fulfill the Poisson equation implying 



<t>k 
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(4) 



Equations ([3]) and <j4j define together a line in a space spanned 
by Qj and the wavelength of the perturbation A = 2n/k ex- 
pressed in units of /l cr it = 47r 2 GErf/K 2 , which separates neutrally 
stable (oj 2 > 0) from exponentially unstable (oj 2 < 0) perturba- 
tions of the disk. The criterion that ensures that all perturbations 
are neutrally stable is the famous Toomre criterion 



6t > 1 



(5) 



The model of a local patch of a galactic disk has been extended 
by Fuchs (2004) by embedding it into a dark halo. All density 
gradients in the halo are neglected as in the disk so that the halo 
density distribution is assumed to be homogeneous. The dark 
matter particles follow straight-line orbits with an isotropic ve- 
locity distribution modelled also by a Gaussian distribution. We 
can directly apply the results of Fuchs (2004). The dark matter 
halo responds to the potential perturbation <l>d in the disk and 
develops potential perturbations <t>h which have the same ra- 
dial structure exp(ikx) as in the disk. From eqns. (26) and (28) 
of Fuchs (2004 ) follows that the Fourier coefficients of the po- 
tential perturbation in the halo at the midplane of the disk are 
given by 



2nGp h 1 



(6) 



where pj, and en, denote the density of the dark halo and the 
velocity dispersion of the dark matter particles, respectively. 

This induced perturbation of the gravitational potential of 
the dark halo has to be taken into account on the rhs of eq. ||3), 



<I>dk + <E> hk °c 4> 



(7) 



which means that the halo supports the perturbation of the disk 
and the density perturbation in the disk is stronger than in an 
isolated disk. Combining eqns. (O to (0 leads to an implicit 
equation that describes the line of neutrally stable (oj = 0) per- 
turbations in the space spanned by Q T and A/A cr \ t as in the case 
of an isolated disk, but now modified by the extra term given in 
eqns. © and (0. This can be cast into dimensionless form as 
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Fig. 1. Separation of stable from unstable perturbations of a 
self-gravitating disk embedded in a live dark halo. Qj denotes 
the usual Toomre stability index and A is the wavelength of the 
perturbation measured in units of A c „t- Unstable perturbations 
are located in the parameter space below the dividing lines. 
Lines are shown for values of the /3 parameter, which describes 
the dynamical responsiveness of the dark halo, f3 — 0, 0.2, 0.25, 
and 0.3, respectively. 



with the parameters a = (3.36/2n) - 0.286 and j3 — 
(2nGp h la 2 l ) ■ (2nGY d ) 2 / k 4 . 

In Fig. Q] we illustrate solutions of eq. ([8]l for various val- 
ues of /3. The case ft = reproduces Toomre's classical ( 1964) 
result. The unstable perturbations (oj 2 < 0) are located in the 
parameter space below the line. Thus for Qj > 1 all perturba- 
tions are neutrally stable (or > 0). This is no longer the case, 
if finite values of (3 are considered. The graphs of the solutions 
shown in Fig. Q] always turn upwards at large wavelengths A. 
Thus at large enough wavelengths all perturbations of the disk 
- halo system become unstable and grow exponentially. This 
behaviour is related to the Jeans collapse of the halo compo- 
nent. Its Jeans length is given by A\ = no-^/Gph or 



A 

^crit 



2/3 



(9) 



In real haloes the Jeans length will be of the order of the size 
of the halo or even larger, because otherwise the haloes would 
have collapsed to smaller sizes. Thus the inferred instability of 
the disk - halo system on large scales seems not to occur in real 
galaxies. As can be seen from Fig. 1 stability on small scales 
can be ensured by Qj indices at thresholds which are slightly 
larger than in isolated disks. 

3. Discussion and Conclusions 

As a first application of the stability criterion derived here we 
test the stability of the Milky Way disk and the surrounding 
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Table 1. Local parameters of the Milky Way 



References 



ITt 38 M Q /pc 2 

cr d 40 km/s 

k V2 • 220 km/s/8.5 kpc 

Ph 0.01 M /pc 3 

cr h 220 km/s/ V2 

/Icrit 4.8 kpc 
39 kpc 



(Holmberg & Flynn 2004) 
(JahreiB & Wielen 1997) 

(flat rotation curve) 
(Bahcall & Soneira [L?80T > 

(isothermal sphere) 

(J3 = 0.0078) 



dark halo in the vicinity of the Sun. The local disk and halo 
parameters listed in Table [TJ imply Qt = 2.8 and /3 = 0.0057, 
respectively. If we include in our estimate the cold interstellar 
gas with a local surface density of 4M /pc 2 (Dame 1993) and 
assume a velocity dispersion of the interstellar gas of cr g = 5 
km/s, which leads to a reduced mass weighted effective veloc- 
ity dispersion of the combined stellar and gaseous disks, the 
parameter values change to Qi - 2.2 and p 1 = 0.0078, respec- 
tively. Equation (0 implies that As — 8 /t cl -i t = 39 kpc. Thus the 
Milky Way disk and halo system seems to be very stable. 

In order to explore in what range the /J-parameter of spiral 
galaxies is to be expected, we consider the model of a Mestel 
disk with the surface density = So R 1 embedded in a singu- 
lar isothermal sphere representing the dark halo with the den- 
sity distribution p/, = poR~ 2 - The rotation curve of the model 
galaxy is given by 



v 2 c (R) = v 2 d (R) + vl(R) 



(10) 



with the disk contribution i>?(7?) = 2nGI.dR = const, and 
the halo contribution = AnGphR 2 = const. (Binney & 

Tremaine 1987). From the radial Jeans equation follows imme- 
diately that the velocity dispersion of the dark matter particles 
is given by 



2 + vh, 



(ID 



because the particles are bound by both the gravitational disk 
and halo potentials. We find then 



P = 



1 R 2 



1 



R 2 v 2 + u 2 



4 0,2 + ^2 4 (U 2 + U 2)3 



which implies the maximal value 



0.037 . 



(12) 



(13) 



This means that in realistic halo models its density cannot be 
increased, on one hand, and the velocity dispersion of the halo 
particles lowered, on the other hand, indiscriminately, because 
the halo model has to stay in radial hydrostatic equilibrium. 
Equation ( fT3l implies /ij = 3 .7 ^t C rit- In order to ensure stabil- 
ity at smaller wave lengths the Toomre stability index must be 
larger than Qj > 1.02. We conclude from this discussion that 
embedded galactic disks are not prone to Jeans instabilities, 
provided their Toomre stability index is a few percent higher 
than Qj — 1 . From a practical point of view the destabilizing 
effect of the surrounding dark halo on the Jeans instability of 
the embedded galactic disks seems to be negligible. 
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